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PROPAGATION OF CHAOS, WASSERSTEIN GRADIENT FLOWS 
AND TORIC KAHLER-EINSTEIN METRICS 

ROBERT J. BERMAN, MAGNUS ONNHEIM 


Abstract. Motivated by a probabilistic approach to Kahler-Einstein met¬ 
rics we consider a general non-equilibrium statistical mechanics model in Eu¬ 
clidean space consisting of the stochastic gradient flow of a given (possibly 
singular) quasi-convex N-particle interaction energy. We show that a deter¬ 
ministic “macroscopic” evolution equation emerges in the large N-limit of many 
particles. This is a strengthening of previous results which required a uniform 
two-sided bound on the Hessian of the interaction energy. The proof uses 
the theory of weak gradient flows on the Wasserstein space. Applied to the 
setting of permanental point processes at “negative temperature” the corre¬ 
sponding limiting evolution equation yields a drift-diffusion equation, coupled 
to the Monge-Ampere operator, whose static solutions correspond to toric 
Kahler-Einstein metrics. This drift-diffusion equation is the gradient flow on 
the Wasserstein space of probability measures of the K-energy functional in 
Kahler geometry and it can be seen as a fully non-linear version of various 
extensively studied dissipative evolution equations and conservations laws, in¬ 
cluding the Keller-Segel equation and Burger’s equation. We also obtain a 
real probabilistic (and tropical) analog of the complex geometric Yau-Tian- 
Donaldson conjecture in this setting. In a companion paper applications to 
singular pair interactions are given. 
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1. Introduction 

The present work is motivated by the probabilistic approach to the construction 
of canonical metrics, or more precisely Kahler-Einstein metrics, on complex alge¬ 
braic varieties introduced in [siini, formulated in terms of certain /3—deformations 
of determinantal (fermionic) point processes. The approach in [SJ I10| uses ideas 
from equilibrium statistical mechanics (Boltzmann-Gibbs measures) and the main 
challenge concerns the existence problem for Kahler-Einstein metrics on a complex 
manifold X with positive Ricci curvature, which is closely related to the seminal 
Yau-Tian-Donaldson conjecture in complex geometry. In this paper, which is one 
in a series, we will be concerned with a dynamic version of the probabilistic ap¬ 
proach in mm- In other words, we are in the realm of non-equilibrium statistical 
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mechanics, where the relaxation to equilibrium is studied. As the general complex 
geometric setting appears to be extremely challenging, due to the severe singu¬ 
larities and non-linearity of the corresponding interaction energies, we will here 
focus on the real analog of the complex setting introduced in [3, taking place in 
R" and which corresponds to the case when AT is a toric complex algebraic vari¬ 
ety. As explained in j9] in this real setting the determinantal (fermionic) processes 
are replaced by permanental (bosonic) processes and convexity plays the role of 
positive Ricci curvature/ plurisubharmonicity (see Section 14.31 for some geometric 
background). 

Our main result lTheorem ll.l[) shows that a deterministic evolution equation on 
the space of all probability measures on R” emerges from the underlying stochastic 
dynamics, which as explained below can be seen as a new “propagation of chaos” 
result. The evolution equation in question is a drift-diffusion equation coupled to 
the fully non-linear real Monge-Ampere operator. It exhibits a phase transition at 
a certain geometrically determined critical parameter. It turns out that in the case 
of the real line (i.e. n = 1) this equation is closely related to various extensively 
studied evolution equations, notably the Keller-Segel equation in chemotaxis [48], 
Burger’s equation [481186] in the theory of non-linear waves and scalar conservation 
laws and the deterministic version of the Kardar-Parisi-Zhang (KPZ) equation de¬ 
scribing surface growth m- In the higher dimensional real case the equation can 
be viewed as a dissipative viscous version of the semi-geostrophic equation appear¬ 
ing in dynamic meteorology (see [5^IT] and references therein). Moreover, closely 
related evolution equations appear in cosmology and in particular in Brenier’s ap¬ 
proach to the Zeldovich model used in the early universe reconstruction problem 
[S51I3H11211122] ■ For the corresponding static problem we establish a real analog 
of the Yau-Tian-Donaldson conjecture (Theorem (T^j) which, in particular, yields a 
probabilistic construction of toric Kahler-Einstein metrics (see Section 14.81 for the 
relations to complex geometry). 

As we were not able to deduce the type of propagation of chaos result we needed 
from previous general results and approaches the main body of the paper establishes 
the appropriate propagation of chaos result, which, to the best of our knowledge, 
is new and hopefully the result, as well as the method of proof, is of indepen¬ 
dent interest. As will be clear below our approach heavily relies on the theory of 
weak gradient flows on the Wasserstein space 7^2 (R") of probability measure 
on R" developed in the seminal work of Ambrosio-Gigli-Savare |2], which provides 
a rigorous framework for the Otto calculus [62| . In particular, as in [2| convexity 
(or more generally A—convexity) plays a prominent role. Our limiting evolution 
equation will appear as the gradient flow on 7^2 (R") of a certain free energy type 
functional F. Interestingly, as observed in |7] the functional F may be identified 
with Mabuchi’s K-energy functional on the space of Kahler metrics, which plays 
a key role in Kahler geometry and whose gradient flow with respect to different 
metrics (the Mabuchi-Donaldson-Semmes metric and Calabi’s gradient metric) are 
the renowned Calabi flow and Kahler-Ricci flow, respectively |28|. The regularity 
and large time properties of the evolution equation appearing here will be studied 
elsewhere [T3||^. 

In the remaining part of the introduction we will state our main results: first, 
a general propagation of chaos result assuming a uniform Lipschitz and convexity 
assumption on the interaction energy and then the application to permanental point 
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processes and toric Kahler-Einstein metrics. In the companion paper m we will 
give a more general formulation of the propagation of chaos result, by relaxing 
some of the assumptions (in particular, this will yield sharp convergence results for 
strongly singular repulsive pair interactions when n = 1). 

1.1. Propagation of chaos and Wasserstein gradient flows. Consider a sys¬ 
tem of N identical particles diffusing on the n—dimensional Euclidean space X := 
R" and interacting by a symmetric energy function E^^\xi,X 2 , ■■■■, xn)- At a fixed 
inverse temperature /3 the distribution of particles at time t is, according to non¬ 
equilibrium statistical mechanics, described by the following system of stochastic 
differential equations (SDEs), under suitable regularity assumptions on : 

(1.1) dx^{t) = --^E^^'>{xi,X 2 ,....,XN)dt+ ^J^dBi{t), 

where Bi denotes N independent Brownian motions on R"; the equation is called 
the (overdamped) Langevin equation in the physics literature. In other words, this 
is the Ito diffusion on R" describing the (downward) gradient flow of the function 
E^^'^ on the configuration space perturbed by a noise term. A classical problem 
in mathematical physics going back to Boltzmann and made precise by Kac |46| is 
to show that, in the many particle limit where iV —>■ oo, a deterministic macroscopic 
evolution emerges from the stochastic microscopic dynamics described byO More 
precisely, denoting by 5^ the empirical measures 

( 1 - 2 ) 

the SDEs [01 define a curve 5M{t) of random measures on X. The problem is to 
show that, if at the initial time t = 0 the random variables Xi are independent 
with identical distribution then at any later time t the empirical measure 6N(t) 
converges in law to a deterministic curve git of measures on R” 

(1.3) lim 6N{t) = 

N—^OO 

In the terminology of Kac [IHI (see also inzi) this means that propagation of chaos 
holds at any time t. It should be stressed that the previous statement admits a pure 
PDE formulation, not involving any stochastic calculus (see Section [HD and it is 
this analytic point of view that we will adopt here. 0 

Of course, if propagation of chaos is to hold then some consistency assumptions 
have to be made on the sequence E^^'> of energy functions as N tends to infinity. 
The standard assumption in the literature ensuring that propagation of chaos does 
hold is that E^^'){xi,X 2 , ■■■■,xn) can be as written as 

(1.4) E^^\xuX2,....,xn)=NE{6n) 

for a hxed functional E on the space of 'P{X) of all probability measures on X, 
where E is assumed to have appropriate regularity properties (to be detailed below). 
This is sometimes called a mean field model. By the results in HE EH El HOI it 
then follows that the limit ptdx) with initial data //o(= podx) is uniquely 


^From a differential geometric point of view the SDEs firTI correspond. under the transformation 
p to the heat flow on of the Witten Laplacian of the “Morse function” E, but we 

will not elaborate on this point here. 
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determined and satisfies an explicit non-linear evolution equation on 'P{X) of the 
following form: 

(1-5) ^ - V • {pth[pt]) 

where we have identified /r(= pdx) with its density p and b[p] is a function on 'P{X) 
taking valued in the space of vector fields on X defined as minus the gradient of 
the differential of E on V{X) 

(1.6) b[p\ = -X{dE\^) 

where the differential dE\^^ at p is identified with a function on X, by standard 
duality (the alternative suggestive notation b[p\ = — is often used in the 

literature). In the kinetic theory literature drift-diffusion equations of the form 
11.61 are often called McKean-Vlasov equations. More generally, the results re¬ 
ferred to above hold in the more general setting where the gradient vector field 
— -^E^^\xi,X 2 , ....,X 7 v) on X appearing in equation II.II is replaced by 6[(5Ar] for a 
given vector field valued function v[p\ on V{X), satisfying appropriate continuity 
properties. 

One of the main aims of the present work is to introduce a new approach to 
the propagating of chaos result 11.31 for the stochastic dynamics 11.11 which exploits 
the gradient structure of the equations in question and which applies under weaker 
assumptions than the previous results referred to above. As indicated above our 
main motivation for weakening the assumptions comes from the applications to 
toric Kahler-Einstein metrics described below. In that case there is a functional 
E{p) on 7^(R") such that 

(1.7) ■^E^^\xi,X 2 , ....,XAr) = E{Sn) + o(l), 

for a sequence of functionals E^^'^ which are uniformly Lipschitz continuous in each 
variable separately, i.e. there is a constant C such that 

(1.8) |V,,aW|<C 

and the error term o(l) tends to zero, as N ^ oo (for Xi uniformly bounded). 
Moreover, E^^'> is A—convex on X^ for some real number A, which, by symmetry, 
means that the (distributional) Hessian are uniformly bounded from below for any 
fixed index i : 

(1.9) > A/, 

where I denotes the identity matrix. This implies that there exists a unique solution 
to the evolution equation 11.51 in the sense of weak gradient flows on the space 
V 2 {X) of all probability measures with finite second moments equipped with the 
Wasserstein metric [2]: 


where Tjg is the free energy type functional corresponding to the macroscopic energy 
E{p) at inverse temperature jd : 

Fp{p) = E{p) + ^H{p), 
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and where H{^) is the Boltzmann entropy of /x (see Section \TJ \ for notation). 

Theorem 1.1. Suppose that is a sequence of symmetric functions on 
satisfying the Main A ssumvtions \l.'T[\1.8\ and \l.tA Then, for any fixed positive time 
t, the empirical measure of the system of SDEs flTTl with independent 

initial data distributed according to po G converges in law, as N ^ oo, to 

the deterministic measure pit evolving by the gradient flow on the Wasserstein space 
of the corresponding free energy functional Ep emanating from /ig. 

It should be stressed that the key point of our approach is that we do not need to 
assume that the drift v[p]{x) defined by formula If .61 has any continuity properties 
with respect to p, or x. Or more precisely, even if the TV— dependent drift vector 
field may very well be smooth for any hxed N, we do not assume that it is 
uniformly bounded in N. This will be crucial in the applications to toric Kahler- 
Einstein metrics below. 

We recall that if the drift has suitable continuity assumptions, then the existence 
of a solution to the drift-diffusion equation 11.51 can be established using fix point 
type arguments m- However, in our case we have, in general, to resort to the weak 
gradient flow solutions provided by the general theory in |2] , where the solution pt 
can be characterized uniquely by a differential inequality called the Evolutionary 
Variational Inequality (EVI). As shown in [2] the corresponding solution pt satisfies 
the drift-diffusion equatioi il.bl in a suitable weak sense (as follows formally from the 
Otto calculus [S2]). 

Our approach is inspired by the approach of Messer-Spohn |5H] concerning the 
static problem for the SDEs 11.11 i.e the study of the Boltzmann-Gibbs measure on 
associated to E^^'> at inverse temperature f : 



Zn 




assuming that the normalization function (the partition function) is finite: 


Z]s[ '■= f e ^ dx < oo 

Jx^ 


From a statistical mechanical point of view this probability measure describes the 
microscopic equilibrium distribution at a fixed inverse temperature /3 and it appears 
as the law of the large time limit (for N fixed) of the empirical measures 5M{t). 
We thus get a uniform approach which applies both to the dynamic and the static 
setting and which in the latter case leads to the following generalization of |58) : 


Theorem 1.2. Suppose that E^^'> is uniformly Lipschitz continuous and convex 
and satisfies the following uniform properness assumption: 




i=l 


for some positive constant C. Then the Boltzmann-Gibbs measures corresponding 
to En are well-defined and 

lim —^^log [ e~^^^ ^dx®^ = inf Ffl(> —cx)) 
iV^oo TV/3 -PalR") 
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Moreover, the corresponding free energy functional Fp on 7^2 (®.”) admits a mini- 
mizer fip and if it is uniguely determined , then the corresponding empirical mea¬ 
sures on K” converge in law as N ^ oo to the deterministic measure fap (which, 
as a conseguence, is log concave). 


It should be stressed that the properness assumption in the previous theorem 
which will appear naturally in the setting of toric Kahler-Einstein metrics below, 
corresponds to properness wrt the Wasserstein metric on 7^(IR"), which is thus 
strictly weaker than demanding properness with respect to the Wasserstein 
metric. But using the convexity assumption on E]\f, we will bypass this difficulty 
using Prekopa’s inequality and Borell’s lemma. 

We briefly point out that the previous theorem is also related to previous work on 
lattice spin models such as the Kac model [42] , as well as lattice models for random 
growth of surfaces [39], where the large limit corresponds to the “thermodynamic 
limit”, where the lattice is approximated by a hnite volume lattice. 


1.1.1. Idea of the proof of Theorem \l.l\ and comparison with previous results. The 
starting point of the proof is the basic fact that the SDEs [TTT] on admit a PDE 
formulation: they correspond to a linear evolution fiN(t) of probability measures 
(or densities) on , given by the corresponding forward Kolmogorov equation 
(also called the Fokker-Planck equation). Given this fact our proof of Theorem ll.il 
proceeds in a variational manner, building on [5]: the rough idea to show that the 
any weak limit curve T{t) of the laws 

rw(i) := iSN)*mit) € W), Y = P2(M”) 


is of the form r(t) := where the curve fit in P 2 (R”) is uniquely determined by 
a “dynamic minimizing property”. To this end we first discretize time, by fixing a 
small time mesh r := tj^i — tj and replace, for any fixed N, the curve P ^{1) with 
its discretized version T^^ftj), defined by a variational Euler scheme (a “minimizing 
movement” in De Giorgi’s terminology) as in |45l[2]. We then establish a discretized 
version of Theorem o saying that if, at a given discrete time tj the following 
convergence holds in the Wasserstein metric 


lim rf = S, 

N^oa ^ 


Vi, > 


then the convergence also holds at the next time step tj+i (using a variational 
argument). In particular, since, by assumption, the convergence above holds at 
the initial time 0 it “propagates” by induction to hold at any later discrete time. 
Finally, we prove Theorem ll.ll bv letting the mesh r tend to zero. This last step uses 
that the very precise error estimates established in [2] , for discretization schemes as 
above, only depend on a uniform lower bound A on the convexity of the interaction 
energies. 

Our proof appears to be to rather different from the probabilistic approaches in 
[5T] (and elsewhere) which are based on a study of non-linear martingales and 
the recent PDE approach in [^ . As pointed out above these approaches require a 
Lipschitz control on the drift vector field and hence a two-sided uniform bound 
on the Hessian of the interaction energy , while we only require a uniform lower 
bound. 

It may also be illuminating to think about the convergence of T^ft) towards 
r(t) as a kind of a stability result for the sequence of weak gradient flows on 
V 2 (Y), associated to the corresponding mean free energies, viewed as functionals 
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on V 2 (Y). This situation is somewhat similar to the stability result for gradient 
flows on 'P 2 {H) in [21 [5], where iJ is a Hilbert space, but the main difference here 
is that the underlying space Y is not a Hilbert space, as opposed to the setting in 
nia, which prevents one from directly applying the error estimates in |2] on the 
space 7^2 (H) itself (this analog will be expanded on in the companion paper |11|1. 

1.1.2. Generalizations. Before continuing with the applications of Theorem 11.11 
(and its static analog) to permanental point process and toric Kahler-Einstein met¬ 
rics we want to stress that the assumptions appearing in Theorem 11.11 may certainly 
be weakened: 

• By rescaling we may as well allow the “inverse temperature” P ap¬ 

pearing in the SDEs 11.11 to depend on N as long as 

Pn —>■ /3 £ [0, oo], 

as TV oo. In particular. Theorem 11.11 also applies to /3 = oo where the 
evolution equation 11.51 becomes a pure transport equation (i.e. with no 
diffusion). However, the precise relation to weak solutions becomes much 
more subtle and is closely related to the notions of entropy solutions and 
viscosity solutions studied in the PDE-literature m (as detailed in my 
In fact, one may even allow that /3jv = oo, where the corresponding con¬ 
vergence results yields a deterministic mean field particle approximation. 

• The assumption 1 1. 71 in conjunction with the Lipschitz assumption 11.81 may 

be replaced by the assumption that the limit of the mean energies cor¬ 
responding to in the sense of statistical mechanics, exists (i.e. that 

Proposition 12 . 201 below holds) and that E^^'> has a uniform coercivity prop¬ 
erty. For example, one can add to the Lipschitz function E^^'^ any term of 
the form NV{Sn), for a given coercive A—convex V functional on P 2 (®.") 
(one then replaces E{iY) with £’(/r) -|- V(/i)). 

• The convexity assumption on E^^'> may be replaced by a generalized con¬ 
vexity property of the corresponding mean energy functional on P 2 (R”^)‘^^ • 

These generalizations will be developed in the companion paper m- 

1.2. Applications to permanental point processes at negative temper¬ 

ature and toric Kahler-Einstein metrics. Let P be a convex body in R" 
containing zero in its interior and denote by Pi the lattice points in P, i.e. the 
intersection of the convex body P with the integer lattice Z". We fix an auxiliary 
ordering pi, ...,pn of the N elements of Pi. Given a configuration (xi,..., xn) of N 
points on X we denote by Per(a:i, ...,xn) the number defined as the permanent of 
the rank N matrix with entries Aij := : 

(1.11) Per(a;i,...,a:Ar) := Per (e^vP^) = , 

where Sn denotes the symmetric group on N letters. This defines a symmetric 
function on which is canonically attached to P (i.e. it is independent of the 
choice of ordering of Pz)0 We will consider the large N limit which appears when 


^In many body quantum mechanics Perfxi,..., xjy) appears as the Af—particle wave function 
for a bosonic system of N particles represented by the N wave functions , i.e. N planar waves 
with imaginary momenta proportional to pj. 
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P is replaced by the sequence kP of scaled convex bodies, for any positive integer 
k. In particular, N depends on k as 


Nk = 


k^ViP) 


n\ 


r(fc-), 


where V{P) denotes the Euclidean volume of P. In this setting the interaction 
energy is defined by 


( 1 . 12 ) 


1 


= ^logPeT{xi,...,XNk) 


To simplify the notation we will often drop the explicit dependence of N on k. 
By the results in [S] the assumptions in Theorem 1 1.1 1 hold with 


Eifi) := 


where C(/r) is the Monge-Kantorovich optimal cost for transporting /r to the uni¬ 
form probability measure vp on the convex body P, with respect to the standard 
symmetric quadratic cost function c{x,p) = —x ■ p. Hence, the corresponding free 
energy functional may be written as 

(1.13) Fp{p) = -C{p) + ^H{p) 


Theorem 1.3. Assume that /3 > 0. Then, for any fixed positive time t, the em¬ 
pirical measure ^ Sx, of the stochastic process 11.11 driven by ll.lill with initial 
independent data distributed aceording to a fj,Q G V 2 QSP) converges in law to the 
deterministic measure pt = ptdx evolving by the gradient flow on the Wasserstein 
space, defined by the functional Fp (formula M.lA) and satisfying the evolution PDE 
in the distributional sense: 

(1-14) ^ + V • {ptVfit) 

where fitix) is the unigue convex function on'ET solving the Monge-Ampere eguation 

(1-15) y^det(a2<^t) = Pt 

(in the weak sense of Alexandrov) normalized so that (j){0) = 0 and satisfying the 
growth condition <f>{x) < 4>p{x), where 4ip{x) := sup^gpp • x. 


Integrating twice reveals that the stationary equation corresponding to the evo¬ 
lution PDE in the previous equation may be written as follows in terms of the 
convex “potential” 0 : 

(1.16) det(92^) = 

where ptdx = pdx := MA^cf)). As shown in (generalizing the seminal result 
in [70]) there is a solution (j) := fip to the previous equation iff the origin is the 
barycenter bp of P, i.e. iff 6p = 0 and then the solution is smooth (see also |20| for 
generalizations). Moreover, the additive group ME acts faithfully by translations on 
the solution space. Note that up to replacing P with fi~^P we may as well assume 
that fi = \ and the corresponding static equation 

(1.17) det(52^) = e-^ 

is precisely the Kabler-Einstein eguation for a toric Kahler potential (f> oi a Kahler- 
Einstein metric with positive Ricci curvature on the toric variety Xp corresponding 
to P, in the case when P is a rational polytope. More precisely, Xp is a toric log 
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Fano variety and (p corresponds to the Kahler potential of a Kahler-Einstein metric 
with conical singularities along the divisor Xp — C*" “at infinity” - the ordinary 
smooth Fano case appears when the polytope P is a reflexive Delzant polytope 
|7Dl [Ml [Sj ; for some background see Section 14.3.11 

Given the relation to Kahler-Einstein metrics it is natural to ask if Theorem 
fOl has a static analog (as in Theorem [ 121 )? However, it follows from symmetry 
considerations involving the action by translations of the additive group M", that 
the corresponding Boltzmann-Gibbs measure 

describing a permanental point processes at negative temperature, is not even well- 
defined, i.e. the partition function diverges! This is a reflection of the fact 
that the static equation 11.161 has a a multitude a solutions (due to the translation 
symmetry) which from a statistical mechanical point of view is a sign of a first order 
phase transition. However, as we will show, this issue can be bypassed though a 
symmetry breaking mechanism where one introduces a “background potential” V (x) 
with appropriate growth at infinity, dictated by P, i.e. 

(1.18) \V{x) — 4>p{x)\<C, 4>p{x) ■.= s\xpp ■ X 

peP 

and replace the interaction energy E^^'‘\xi, xni,) defined by formula | 1 . 12 l bv 
the convex combination 

(1.19) ^ 

...,xn^) := 7 ilogPer(xi, ...,xn^) + (1 - 7) (^(a^i) + ''' + V{xn^)) 

for a given parameter 7 G [ 0 , 1 ] (which from the point of view of permanental point 
processes plays the role of minus the inverse temperature). Then the convergence 
in Theorem 11.31 still holds with F replaced by 

= -lC{p) -f (1 - 7 ) J Vdp + H{p) 

and the corresponding static equation now becomes 

( 1 . 20 ) detid^P) = 

which has at most one solution for any given 7 G [0,1] and convex body P. Moreover, 
if P satisfies the barycenter condition bp = 0 then there exists a solution to the 
eauation ll. 20 l for any 7 G [ 0 , 1 ] and as 7 —)• 1 it follows from the results in (THIIH] that 
the solutions converges to a particular solution to the Kahler-Einstein equation 
11.171 singled out by the potential V. 


Theorem 1.4. Let P be a convex body in R" containing 0 in its interior and denote 
by bp the barycenter of P. For any potential V in R" satisfying the growth condition 
we have 

• Ifbp = 0, then, for any 7 G [0,1[ the Gibbs measure corresponding to 
the energy function E^fp\xi, ...,X]si,,), he. 


(AT) 

^ 00.7 - 






{Per{xi, ...,xn)) (e '^'>^dx)®^'‘ 


is well-defined and equal to the weak limit, as t ^ 00 of the law of the 
empirical measures for the corresponding SDEsMoreover, as N ^ 00 
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the corresponding empirical measures converge in law to the deterministic 
measure defined as = MA((f>^) for the unique (mod R) solution cj)^ 
of the equation \1.2(A 

• More generally, the Gibbs measure above is well-defined precisely for 7 < 
Rp, where Rp € [0,1] is the following invariant of P : 


( 1 . 21 ) 


Rp := 


l|g|| 


where q is the point in dP where the line segment starting at bp and passing 
through 0 meets dP. Moreover, for any such parameter 7 the corresponding 
convergence results still hold. 


As indicated in the introduction this result can be viewed as a probabilistic 
analog of the seminal Yau-Tian-Donaldson conjecture saying that a Fano manifold 
X admits a Kahler-Einstein metric if and only if X is K-stable (the conjecture has 
very recently been settled by Chen-Donaldson-Sun [261I27| : see also Tian [Ml)- The 
latter notion is of algebro-geometric nature, but in the toric setting it is equivalent 
to the corresponding polytope P having zero as its barycenter (see Section [YS] and 
Corollary 14.II for a comparison with the complex geometric setting). 

From a statistical mechanical point of view the critical value Rp appearing above 
can be seen as a real analog of the well-known critical value appearing in the study 
of the Keller-Segel equation as well as in the study of the 2D log gas in [251149| . 
This connection will be further expanded on elsewhere na, but the main point is 
that the invariant Rp may also be characterized as the sup over all 7 €]0, oo[ such 
that the free energy type functional is bounded from below (compare |S]). 

As we point out in Sections l3.51l4.3?T] our results also apply to the tropical analog 
of the permanental setting above, which can be viewed as the tropicalization of the 
complex geometric setting on the corresponding toric variety. In the corresponding 
deterministic setting (i.e. Pn = 00 ) the particles then perform zigzag paths in R” 
generalizing the extensively studied Sticky Particle System on M [631 [2^ iBl] . This 
is closely related to the Zeldovich model for the formation of large-scale structures 
in cosmology; see [381 [?I1 \T2\ (compare the discussion in Section ing . 


1.3. Acknowledgments. It is a pleasure to thank Eric Carlen for several stim¬ 
ulating discussions and whose inspiring lecture in the Kinetic Theory seminar at 
Chalmers concerning m drew our attention to the Otto calculus and the theory of 
gradient flows on the Wasserstein space. Thanks also to Luigi Ambrosio for helpful 
comments on the paper and to Yann Brenier, Jose Carrillo, Maxime Hauray and 
Bernt Wennberg for providing us with references. This work was supported by 
grants from the Swedish Research Council, the Knut and Alice Wallenberg Foun¬ 
dation and the European Research Council. 

1.4. Organization. In Section we start by recalling the general setup that we 
will need from probability, the theory of Wasserstein spaces and weak gradient 
flows and then turn to the proof of Theorem 11.11 in Section 12.61 (starting with the 
discretized situation). Then the proof of the corresponding static result. Theorem 
11.21 is given. In Section |3l we go on to apply the previous general results to the 
permanental setting, as introduced in Section [3] and its tropical analog. In the 
final section we provide on outlook on some relations to conservation laws, sticky 
particle type systems and complex geometry. The appendix recalls the basics of the 
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formal Otto and is included to serve as a motivation for the material on Wasserstein 
gradient flows. The rather lengthy setup and preparatory material in Section [5] is 
due to our effort to make the paper readable to a rather general audience. 


2. General setup and proof of Theorem 11.11 

2.1. Notation. Given a topological (Polish) space Y we will denote the integration 
pairing between measures ^ on Y (always assumed to be Borel measures) and 
bounded continuous functions / by 


{/,m) := J /m 

(we will avoid the use of the symbol d/r since d will usually refer to a distance 
function on Y). In case Y = we will say that a measure fi has a density, 
denoted by p, if p, is absolutely continuous wrt Lebesgue measure dx and p = pdx. 
We will denote by the space of all probability measures and by Vaci^^) 

the subspace containing those with a density. The Boltzmann entropy H{p) and 
Fisher information I{p) (taking values in ] — oo,oo]) are defined by 

(2.1) H{p) := f {logp)pdx, I{p)= [ ^-^^dx 

Jrx> Jrd P 

(assuming that Vp G L^{dx) and p~^'Vp G Lf{pdx)). More generally, given a 
reference measure pQ on Y the entropy of a measure p relative to po is defined by 



if the probability measure p on W is absolutely continuous with respect to p and 
otherwise F[{p) := oo. The relative Fisher information is defined similarly. Given a 
lower semi-continuous {Isc, for short) function V on Y and (3 g]0, oo] (the “inverse 
temperature) we will denote by Fj the corresponding (Gibbs) free energy functional 
with potential V : 

(2.3) F^{p)-.= I^Vp+^H^,ip), 

which coincides with ^ times the entropy of p relative to e~^po. In particular, since 
Hfj.o{p) > 0, when p and po are probability measures, with equality iS p = po (by 
Jensen’s inequality) the Boltzmann-Gibbs measure 


P-/3V 


Z := 


/ 


-V 

e Po 


of y, at inverse temperature (3 is the unique minimizer of Fj on the space VfY) 
of probability measures, under the integrability assumption that Z < oo (this is 
usually called Gibbs’ variational principle). 


2.2. Wasserstein spaces and metrics. We start with the following very general 
setup. Let {X, d) be a given metric space, which is Polish, i.e. separable and 
complete and denote by 'P{X) the space of all probability measures on X endowed 
with the weak topology, i.e. pj —>■ p weakly in V{X) iff pjf —>■ p/ for any 
bounded continuous function f on X (this is also called the narrow topology in the 
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probability literature). The metric d on X induces type metrics on the iV—fold 
product for any given p S [ 1 , oo[: 

N 

dp{xi, XN] yi, ■■■, Vn) ■■= d{xi,yi)py/p 

The permutation group Sn on TV—letters has a standard action on , defined 
by (cr, (a;i, ccat)) {xa-{i), ■■■,Xa-{N)) and we will denote by X^^'> and tt the 
corresponding quotient and quotient projection, respectively: 

(2.4) := X^/S^, tt: X^ ^ X^^^ 

The quotient X^^'> may be naturally identified with the space of all configurations 
of N points on X. We will denote by d(p) the induced distance function on X^^\ 
suitably normalized: 


1 

dx(.N)^ip{xi,...,XN;yi,.-,yN) ■= 

The normalization factor I/N^^p ensures that the standard embedding of X^^'> into 
the space 'P{X) of all probability measures on X : 

(2.5) xW-^iP(X), (xi,..,cr^)^d^:=1^4. 

(where we will call 5n the empirical measure) is an isometry when V{X) is equipped 
with the LP— Wasserstein metric dwp induced by d (for simplicity we will also write 
dwp = dp) : 


(2.6) d^ {fj.,v)--ini f d{x,y)P-f, 

'f JxxX 

where 7 ranges over all couplings between y, and v, i.e. 7 is a probability measure 
on X X X whose first and second marginals are equal to y and v, respectively (see 
Lemma [2.31 belowl. We will denote W^(X, d) the corresponding LP—Wasserstein 
space, i.e. the subspace of V{X) consisting of all y with finite p th moments: for 
some (and hence any) Xq G X 

/ d{x, xo)Py < 00 

Jx 

We will also write WP{X,d) = 'Pp{X) when it is clear from the context which 
distance d on X is used. 


Remark 2.1. In the terms of the Monge-Kantorovich theory of optimal transport 
|69| d^ {y, v) is the optimal cost to for transporting y to v with respect to the 
cost functional c{x,p) := d{x, y)P. Accordingly a coupling 7 as above is often called 
a transport plan between y and iz and it said to be defined be a transport map 
T if 7 = (/ X r),p where T*/r = v. In particular, if X = R”, p = 2 and y and 
v are absolutely continuous with respect to Lebesgue measure, then, by Brenier’s 
theorem nni, the optimal transport plan 7 is always defined by a transport map 
T(:= Tp ) of the form = V^, where ^ is a convex function on K" (optimizing 
the dual Kantorovich functional). 

We recall the following standard 
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Proposition 2.2. A sequence nj converges to n in the distance topology in W^{X, d) 
iff fXj converges to p, in the weakly in V{X) and the p th moments converge (the 
latter assumption is automatic if X is compact). In particular, if fij converges to 
fi weakly in V{X) and the p th moments are uniformly bounded, then Hj converges 
to fi in the distance topology in {X, d) for any p' < p. 

Proof. For the first statement see for example [691 Theorem 7.12]. The second 
statement is certainly also well-known, but for completeness we include a simple 
proof. Decompose 


/ d{x,xoy = / 
Jx J{d 


d{x,xoy fJ-j + / d{x,xoy fJ-j 

f {cZ(a:,3;o):<^} {d(a:,3;o)>^} 

By the assumption and Chebishev’s inequality the second terms may be estimated 
from above by C/R^p~p'^ and by the assumption of weak convergence the first 
term converges to d{x,xoY fi as j ^ oo (by taking / to be a suitable 

regularization of I{d(x,xo)<R}dix, xqY )■ Finally, letting R tend to infinity concludes 
the proof. □ 

Since Yp := (Wp{X),dwp)i-= 'Pp(X)) is also a Polish space we can iterate the 
previous construction and consider the Wasserstein space Wq{Y) C V{V{X)) that 
we will write as Wq{Vp{X)), which is thus the space of all probability measures F 
on V{X) such that, for some po G Wp{X) 


L 


dp{pL,pLoyT < oo 


V(X) 

Lemma 2.3. (Three isometries) 

• The empirical measure Sn defines an isometric embedding {X^^\ d(^py 

T^piX) 

• The corresponding push-forward map {5 n)* from V{X^^) to V(V{X)) in¬ 
duces an isometric embedding between the corresponding Wasserstein spaces 
W',(XW,d(p)) andWqiVpiX)). 

• The push-forward tt* of the quotient projection tt : X^ —> X^^'> induces an 
isometry between the subspace of symmetric measures in {Wq{X^, j^rj^dp) 
and the space {Wq{X^^\ d(^p-j) 

Proof. The first statement is a well-known consequence of the Birkhoff-Von Neu¬ 
mann theorem which gives that for any symmetric function c{x, y) on XxX we have 


that if p = i Y.i=i dxi and v = Ya=i dyi forgiven (xi, ...,xn), (yi,..., i / at ) € X 
then 


N 


inf / c{x,y)dY = inf / c{x,y)dT 

J ^n(u,v)J 


where F v) C r(p, v) consists of couplings of the form Fo- X) dxi ® , 

for a € ^jv, where S'jv is the symmetric group on N letters. The second statement 
then follows from the following general fact: if / : (Yi,di) —>■ ( 12 ,^ 2 ) is an isome¬ 
try between two metric spaces, then /* gives an isometry between Wq(Yi,di) and 
Wg(T 2 , df). This follows immediately from the definitions once one observes that one 
may assume that the coupling 72 between /»p and is of the form /* 7 i for some 
coupling 71 between p and v. The point is that 7 can be taken to be concentrated 
on /(hi) X /(F 2 ) (since this set contains the product of the supports of p and v) and 
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hence one can take 71 := (/ ^ 0 / ^)*72 where (/ ^ ® / ^)if{y)Jiy')) ■= {y,y') 
is well-defined, since / induces a bijection between Yi and /(hi). Finally, the last 
statement follows immediately from the following general claim applied to h" = 
with d = j^l/p dx^ ip and G = Sx- Let G be a compact group acting by isometries 
on a metric space {Y,d) and consider the natural projection ir : Y ^ Y/G. We 
denote by da the induced quotient metric on Y/G. The push-forward tt* gives a 
bijection between the space V{X)^ or all G—invariant probability measures on X 
and 'P{X/G). The claim is that tt* induces an isometry between the corresponding 
Wasserstein spaces Vq{X)^ and Vq{X/G) i.e. dwq{y-,v) = du/^(7r*/r,7r*i/) if and 
v are G—invariant (see |581 Lemma 5.36] Lemma 5.36). □ 

Let us also recall the following classical result 

Lemma 2.4. Let /tq be a probability measure on X. Then —>■ in 

'P{V{X)) weakly as TV 00 

In fact, according to Sanov’s classical theorem the previous convergence results 
even holds in the sense of large deviations at speed N with rate functional given 
by the relative entropy functional H^g{-) [321 Theorem 6.2.10]. We note that that 
a (somewhat non-standard) proof of this classical result can be obtained using the 
argument in the proof of Theorem 11.21 applied to E = 0. 

2.2.1. The present setting. We will apply the previous setup to X = R" endowed 
with the Euclidean metric d. Moreover, we will mainly use the case p = 2. Then the 
corresponding metric d 2 on X^ is the Euclidean metric on X^ = Identifying 
a symmetric (i.e. invariant) probability measure y-x on X^ with a probability 
measures on the quotient X^^'> (as in Lemma (2.31) the second and third point in 
Lemma 12.31 may (with q = 2) he summarized by the following chain of equalities 
that will be used repeatedly below: 

(2.7) —d2{yN, y'xY = d(2){yN, = c? 2 (rAr,r(y)^, 

where T x and r(y denote the push-forwards under 5x of yx and , respectively. 

2.3. The Main Assumptions on the interaction energy Set X = R” 

and denote by d the Euclidean distance function on X. Throughout the paper E^^'> 
will denote a symmetric, i.e. S'jv—invariant, sequence of functions on X^ and we 
will make the following Main Assumptions: 

(1) The functional E^^'> is uniformly Lipschitz continuous {X,d) in each vari¬ 
able, or equivalently, under the isometric embedding of X^^'> in 'P{X) (by 
the empirical measure 5x) the sequence E^^'>/N extends to define a se¬ 
quence of functionals on 7 ^ 2 (A) which are uniformly Lipschitz continuous. 

(2) The sequence E^^'> /N of functions on 7^2 (A) has a unique point-wise limit 
E{y). 

(3) The sequence E^^'> is A—convex on (A, d), uniformly in N. 

Note that since E{y) above is assumed Lipschitz continuous on 7 ^ 2 (A) it is uniquely 
determined by its restriction to the subspace V{X)c consisting of all y with compact 
support. 

Lemma 2.5. Assume given a sequence E^^'> satisfying the uniform Lipschitz as¬ 
sumption (1). 
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• Then the second point is (2) is equivalent to point-wise convergence of 
E^^')/N towards E{p) for any /x in V{X) with compact support 

• The second point (2) implies that 

(2.8) lim 1 = E{^) 

N^oo I\ J 

for any pL e V 2 {X). 

• Conversely, if \2.8\ holds for any /x G 'PiX) with compact support then 
E^^'>/N converges towards E{p,) for any /x in V 2 {X). 

Proof. Given p, G V 2 iX) we define the truncation pn := ^ gy gjp, 

J ^{d(x,XQ)<RM 

assumption /N)(p) — (E^^^/N)(pji)l < Cd 2 {p, pr). In particular, on '■= 

{E^^'> /N){p) is a uniformly bounded sequence in R. Next, letting N ^ oo gives 
\a — E{pr)\ < d2{p, Pr) for any limit point a G M of the sequence gn- Finally 
letting R ^ oo and using the Lip continuity forces a = E{p), as desired and hence 
{E^^'>/N){p) converges towards E{p), as desired. To prove formula |2.8I we first 
remark that it follows from the general convergence in Proposition 12. 201 below that 

hm 1 f E<-^^p^^ = f E{,2)ri,2), 

N-fco N J Jv{X) 

where F is a weak limit point of {5m)*P®^ ■ But by Lemma 12.41) the limit point is 
unique and given by F = 6^. Hence, the rhs above is equal to E(p), as desired. □ 

2.4. The forward Kolmogorov equation for the SDEs and the mean free 
energy Fn. Fix a positive integer N and /3 > 0 (which may depend on N when 
we will later on let N —>■ oo). Let {X,g) be a Riemannian manifold and denote by 
dV the volume form defined by g. In our case {X, g) will be the Euclidean space 
R". As is well-known, under suitable regularity assumptions, the SDEs ITT] on X^ 
define, for any fixed T, a probability measure tjt on the space of all continuous 
curves (“sample paths”) in X^, i.e. the space of continuous maps [0,T] —>■ X^ (see 
for example inn and reference therein). For t hxed we can thus view x^^^t) as a 
valued random variable on the latter probability space. Then its law 

p[^^ ■= {x^^\t))^r]t 

gives a curve of probability measures on X^^'' of the form p[^'^ = p^^^dV, where 
the density p\^'^ satisfies the corresponding forward Kolmogorov equation: 

(2.9) ^ + V • (p^VEW), 

which thus coincides with the linear Fokker-Planck equation 15.61 on X^ with po¬ 
tential V := E^^\ In particular, the law of the empirical measures SN{t) for the 
SDEs ITT] can be written as the following probability measure on V{X^) : 

Tj^ft) := {5n)*p[^\ 

where Sn is the empirical measure defined by formula 12.51 

Anyway, for our purposes we may as well forget about the SDEs 11.11 and take 
the forward Kolmogorov equation 12.91 on X^ as our the starting point. We will 
exploit the well-known fact, going back to [45] (see Prop |2T4] below) that the latter 
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evolution equation can be interpreted as the gradient-flow on the Wasserstein space 
W 2 {Y), of the functional 


Jx^ P 

where H{-) is the entropy relative to /tq := 1formula 12.21) : occasionally we 

will omit the subscript j3 in the notation 

Following standard terminology in statistical mechanics we will call the scaled 
functional Fm := F^^'>/N the mean free energy, which is thus a sum of the mean 
energy Fm{pn) and the mean entropy F[m{P'n) '■ 


Fm = Fm + sHm, 


(2.10) Fm{pm) ■= j;^F^^^p,m) = ^ 

Note that it follows immediately from the definition that the mean entropy is ad¬ 
ditive: for any p. € V^X) 


= Hit,) 

In case dV is a probability measure it follows immediately from Jensen’s inequal¬ 
ity that iJ(p) > 0. In our Euclidean setting this is not the case but using that 
/ I dx < oo for any given e > 0 one then gets 

( 2 . 11 ) H{p)>-eJ\x\‘^t,-C, 

As a consequence we have the following 

Lemma 2.6. If the mean energy satisfies the uniform coercivity property 

(2.12) 1/ F;W(p^)>-^d2(Miv,r,)^-C 

tv J 2 t* 

for some fixed r, > 0 and F, G W 2 {V{X)) and positive constant C, then so does 
fW/at. 

For example, this is trivially the case under a uniform Lipschitz assumption on 
FW 

(as in the Main Assumptions). 

Remark 2 . 7 . The linear forward Kolmogorov equation 12.91 can also be viewed as 
the gradient flow of the mean free energy jf Qj^g instead uses the scaled 

metric gx '■= on X^. Moreover, in our case F^^v) i^g symmetric, i.e. 

Fiv—invariant and hence the flow defined wrt {X^,gM) descends to the flow defined 
with respect to /Sm equipped with the distance function dx(N) defined 

in Section 12.21 Using the isometric embedding defined by the empirical measure 
(Lemma 12 . 31 ) we can thus view the sequence of flows on the sequence of spaces 
V{X^) as a sequence of flows on the same (infinite dimensional) space W 2 {V{X)) 
and this is the geometric motivation for the proof of Theorem 11.11 
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2.5. Gradient flows on the L^—Wasserstein space and variational dis¬ 
cretizations. In this section we will recall the fundamental results froni[2] that 
we will rely on. Let G be a lower semi-continuous function on a complete metric 
space (M, d). In this generality there are, as explained in [5], various notions of weak 
gradient flows ut for G (or “steepest descents”) emanating from an initial point uq 
in M, symbolically written as 

(2.13) '^=-' s / G { ut ), limu(t) = Mo 

dt t-s-o 

The strongest form of weak gradient flows on metric spaces discussed in [5] con¬ 
cern A—convex functionals G and are defined by the property that ut satisfies the 
following Evolution Variational Inequalities (EVI) 

(2.14) 

^^d^{ut,v)+G{u{t))+ ^d^{^jbt,vY <G{v) a.e. t > 0, Vv £ M : G{v) < oo 

together with the initial condition lim(_>.ou(t) = uq in (M,d). Then ut is uniquely 
determined by ug, as shown in [21 Cor 4.3.3] and we shall say that ut is the EVI- 
gradient flow of G emanating from ug. We recall that A— convexity on a metric space 
essentially means that the distributional second derivatives are bounded from below 
by A along any geodesic segment in M (compare below). When M has Non-Positive 
Curvature, NPC (in the sense of Alexandrov) the existence of a solution ut satisfying 
the EVI was shown by Mayer EH for any lower-semicontinuos A—convex functional, 
by mimicking the Crandall-Liggett technique in the Hilbert space setting. 

However, in our case {M,d) will be the Wasserstein space P 2 (R^) for the 
space of all probability measures /r on which does not have non-positive cur¬ 
vature (when d > 1). Still, as shown in |2], the analog of Meyer’s result does hold 
under the stronger assumption that G be A—convex along any generalized geodesic 
/is in For our purposes it will be enough to consider Isc A—convex func¬ 

tionals G with the property that 7^2,ac(R‘^) is weakly dense in {G < oo}. Then the 
A—convexity of G means (compare [2] Proposition 9.210]) that for any generalized 
geodesic /ig = Psdx in P 2 ,ac(R‘^) the function G{ps) is continuous on [0,1] and the 
distributional second derivatives on ]0,1] satisfy 

d^Gjp.) ^ , 
d^s - ’ 

We recall that a generalized geodesic ps connecting pg and pi in P 2 ,ac(R‘^) is de¬ 
termined by specifying a “base measure” i/ G 'P 2 ,ac(R‘^))- Then ps is defined as the 
following family of push-forwards: 

Ps = ((1 - s)Tg + sTi)^ V 

where Ti is the optimal transport map (defined with respect to the cost function 
]a: — 2 /P/ 2 ) pushing forward v to pi (compare Remark l2.1|l . 

Remark 2.8. The bona fide Wasserstein geodesics in P 2 ,oc(K‘^) are obtained by 
taking v = pg (the study of convexity along such geodesics was introduced by 
McCann [54], who called it displacement convexity). But as shown in |2] the point 
of working with general base measures v is that they can be adapted to the discrete 
variational scheme for constructing EVI-gradient flows by taking v = pt^ at the jth 
time step (compare Section [2311]) . 
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We will be relying on the following version of Theorem 4.0.4 and Theorem 11.2.1 
in 12 ]: 


Theorem 2.9. Suppose that G is a Isc real-valued functional on 7^2(R^) which is 
X—convex along generalized geodesics and satisfies the following coercivity property: 
there exist constants r*,C > 0 and pL^ G 7^2(R'^) such that 

(2.15) G{-)>--d2{-,pi,f-C 

r* 

Then there is a unique solution pt to the EVTgradient flow of G, emanating from 
any given /tq G {G < oo}. The flow has the following regularizing effect: pt G 
{|9G| < oo} C {G < oo}. Moreover, G{pit) o,nd e'’'* |(?G| (pt) o-re decreasing, where 
|c7G| denotes the metric slope of G : 


|aG| in) ■- limsup 

iy—¥fL 


{G{:z) - G{p)r 
d(/i, v) 


Remark 2.10. Many more properties of the EVTgradient flow pt are established in 
|2|. For example, pt defines an absolutely continuous curve R ^ 7^2 (K") (in the 
sense of metric spaces) which is locally Lipschitz continuous on ]0, oo[ which is a 
A—contracting semi-group. Moreover, the flows are stable under suitable approxi¬ 
mation of the initial data and the functional G. 


Under suitably regularity assumptions it shown in [5] that the EVTgradient flow 
/it = ptdx furnished by the previous theorem satisfies Otto’s evolution equation 
(recalled in the appendix) in the weak sense: 

Proposition 2.11. Suppose in addition to the assumptions in the previous theorem 
that pLt has a density pt for t > 0. Then pt satisfies the the continuity equation ] 5. 5\ 
in the sense of distributions on x R with 

vt = -{d°G){ptdx), 

where d^G denotes the minimal subdifferential of G. 

We recall that under the assumptions in the previous theorem (and assuming 
{|9G|^ < oo} C 7 ^ 2 ,ac(R”)) the many-valued subdifferential dG on the subspace 
7 ^ 2 ,ac(R") is a metric generalization of the (Frechet) subdifferential Hilbert space 
theory; by definition, it satisfies a “slope inequality along geodesics”: 

( 5 G)(/i) := G L^ip) : 'iv : G{v) > G(/x) -h (C,^')L2(^) + ^d2{v,pf, v{x) := Tf^{x) -x 

where Tff denotes the optimal transport map between p, and iz, as in Remark 12.11 
(note that v is the tangent vector field at 0 of the geodesic ps from p to v). The 
minimal subdifferential d^G on 7 ^ 2 ,ac(®.”) at p is defined as the unique element in 
the subdifferential dG at p minimizing the norm in Lf{p)] in fact, its norm 
coincides with the metric slope of G at /t (in [2| there is also a more general notion 
of extended subdifferential which, however, will not be needed for our purposes). 

Example 2.12. In the case when G = 77 is the Boltzmann entropy and p satisfies 
i7(/i) < oo, so that p has a density p, we have {d°H){p) = p~^S/p G L'^{p) and 
hence 

\dH\^{p)=Iip) 

is the Fisher information of p ('formula l2.ip : see [U Theorem 10.4.17] 
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The following result goes back to McCann (see also [5] for various elabora¬ 
tions) : 

Lemma 2.13. The following functionals are Isc and X—convex along any general¬ 
ized geodesics in 7^2 

• The “potential energy” functional V(/r) := fVfi, defined by a given Isc 
X—convex and Isc function V on (and the converse also holds) 

• The functional pi ^ ^ defined by a given X—convex function Vn on 

and in particular the “interaction energy” functional 

W{p) := j W{x - y)pL{x) ® pl{x) 

defined by a given Isc X—convex function W on 

• The Boltzmann entropy H(fi) (relative to dx) is Isc and convex along any 
generalized geodesics. 

In particular, for any X—convex function V on the corresponding free energy 
functional (formula \2.S\) is X—convex along generalized geodesics, if (3 S]0, oo]. 

Combining the results above we arrive at the following 

Theorem 2.14. Assume given (3 g] 0, oo]. Let E{fi) he a Isc functional on 7^2(R‘^) 
which is X—convex along generalized geodesics and satisfies the coercivity condition 
\2.15\ . Then the EVI-gradient flow pit of the corresponding free energy functional 
Fp := E F H/ (3 exists. Moreover, if (3 < oo, then pit = Ptdx, where pt has finite 
Boltzmann entropy. In particular, 

• If V is a Isc finite X—convex function on M.‘^, then the gradient flow of F^ 
exists (defining a weak solution of the corresponding forward Kolmogorov 
equation/Fokker-Planck equation) 

• If E{p) is a Lipschitz continuous functional on which is X—convex 

along generalized geodesics, then the gradient flow exists for any initial data 
P-o G 7^2(R") o,nd if (3 < oo, then pit = ptdx, where pt has finite Boltzmann 
entropy and Fisher information and the following continuity equation holds 
in the distributional sense on R" x R 

(2.16) ^ -k V(ptut), 

where Vt = d^E is the minimal subdifferential of E at pt = ptdx. 

Proof. By the previous Lemma Fp is also Isc and A—convex and by Lemma |2.6I it 
also satisfies the coercivity condition. Hence, the EVI-gradient flow exists according 
to Theorem 12.91 Moreover, by the general results in [2] Fp is decreasing along the 
flow and in particular locally uniformly bounded from above on ]0, oo[. But, by the 
coercivity assumption E > — oo on P 2 (®.‘^) and hence it follows that II{pt) < oo. 
The second statement then follows by the previous lemma and the fact that the 
coercivity condition holds: by A—convexity f{x) := r;(a;)-|-A|a:p is convex and hence 
fix) > —C'|a;| for some constant C, proving coercivity of v. To prove the last point 
first observe that E{p) > —A — Bdip, p/fi < oo on P 2 (R") by the Lip assumption. 
Since Fpipt) < C it follows that Hipt) < oo, which in particular implies that pt 
has a density pt. Moreover, by Theorem 12.91 \dFpipt)\ < oo for t > 0. But since 
E is assumed Lip continuous we have \dEpipt)\ < oo iff \dHipt)\ < oo, which 
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means that I{^t) has finite Fisher information (see Example i l2.12l) . Finally, the 
distributional equation follows from Proposition l2.11l □ 


2.5.1. The variational discretization scheme (“minimizing movements”). We recall 
that the proof of Theorem 12.91 in |2] uses a discrete approximation scheme intro¬ 
duced by De Giorgi, called the minimizing movement scheme. It can be seen as a 
variational formulation of the (back-ward) Euler scheme. Consider the fixed time 
interval [0,T] and fix a (small) positive number r (the “time step”). In order to de¬ 
fine the “discrete flow” uj corresponding to the sequence of discrete times tj := jr , 
where tj < T with initial data uq one proceeds by iteration: given Uj € M the next 
step Uj+i is obtained by minimizing the following functional on {M, d) := W 2 (K.'^) : 

u I—)• —d{u, UjY -\- G{u) 

2t 

Next, one defines u^{t) for any t G [0,T] by setting u^(tj) = mJ and demanding 
that {t) be constant on ]tj,tjj.i[ and right continuous (we are using a slightly 
different notation than the one in Chapter 2]). 

The curve Ut is then defined as the large m limit of in (M, d) ; as shown in [5] 
the limit indeed exists and satisfies the EVI 12.141 and is thus uniquely determined. 
More precisely, the following quantitative convergence result holds (see [21 Theorem 
4.07, formula 4.024] and [21 Theorem 4.09]): 

Theorem 2.15. Let G be a functional on 7^2(R") satisfying the assumptions in 
Theorem, \2. ,91 with A > 0. Then 

d^{u^{t),u{t)) < i]rnaGp(uo), 

where ]9G|(uo) denotes the metric slope of G at uq. If G is only assumed to be 
X—convex for some, possibly negative, A then 

d{u'^{t), u{t)) < G|r|(G(uo) — inf G), 

for some constant C only depending on A and T. 

Remark 2.16. By the last paragraph on page 79 in [2] even if A < 0 one does not 
need a lower bound on inf G if one replaces |t| with as long as uq is assumed 

to satisfy G{uo) < oo. 

2.6. Proof of propagation of chaos in the discretized setting. In this section 
we fix once and for all the time interval [0,T] and the time step r > 0. We denote 
by the corresponding discretized minimizing movement of the free energy 

functional on 'P 2 {X^,d 2 ) with given initial data The sequence 

is well-defined according to Theorem 12.91 and the Main Assumptions. Moreover, 
by the third isometry property in Lemma 12.31 may be identified with the 
minimizing movement of the mean free energy functional /N on V 2 {X^, d{ 2 )), 
which in turn embeds isometrically to give a discrete flow F^^^ in W 2 {V 2 {X),d 2 ). 
Similarly, we denote by the discretized minimizing movement of the functional 
F on V 2 {X,) with given initial data po- 

Theorem 2.17. Assume that at time tj 

lim (4Ar)*p^f ^ = Sf, 
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in the L'^—Wasserstein metric. Then, at the next time step 

lim ,, 

in the L'^ — Wasserstein metric. 

We recall that given the next measure is defined as the minimizer of 
the following functional on : 

(2-U) Pitk) ■■= 

2.6.1. Proof of Theorem \2.1'A We start with the following direct consequence of 
Proposition 12.21 combined with Lemma 12.31 

Lemma 2.18. Let fXN be a sequence of symmetric probability measures on and 
denote byTjsi := {5 ]si)*P-n the corresponding probability measures on'P{X). Assume 
that the dq —distance o/Pjv to a fixed element in the Wasserstein space Wq{P2{X)) 
is uniformly bounded from above, for some fixed q £ [l,oo[. Then, after perhaps 
passing to a subsequence, there is a probability measure T in Wq{T’ 2 iX)) such that 

lim {6N)*tiN = r 
N ^oo 

weakly in 'P{X) or more precisely in Wq'{'P 2 {X) if 1 < q' < q 

We next recall the following well-known result about the asymptotics of the mean 
entropy (proved in [M|; see also Theorem 5.5 in [41] for generalizations). The proof 
is based on the sub-additivity properties of the entropy. 

Proposition 2.19. Let be a sequence of probability measures on X^ such 

that (5Ar)*/rjv converges weakly <o L £ P(fP{X)). Then 


[ H{p)T 

N^oo Jv(X) 

We will also use the following result, which generalizes a result in [58] concerning 
the case when En is quadratic: 

Proposition 2.20. Let be a sequence of probability measures on X^ such 

that Pjv := {Sn)*Pn converges to P in Wi{V 2 [X)). Then 

hm 1 f = [ Eip)T 

N^oo J d'P(x) 

Proof. Recall that the Wasserstein distance di on P(T), where Y = 'P 2 {X), 
admits the following dual representation: 

d{p, v) = sup / u{p — v) 
u^Lipi J 

where u ranges over all Lip-functions on Y with Lip-constant one. By assumption 
(2.18) di(rjv,r)^o. 
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Using the empirical measure 6n we identify with a uniformly Lipschitz 

continuous sequence of functions on V{X) which by the Main Assumptions point- 
wise to to E{fj.). First observe that since is uniformly Lipschitz continuous 

we have that 

lim f N-^E’-^^Tn -T) =0 

using I^TSl combined with the dual representation of the L^—Wasserstein distance. 
Hence, 

lim — [ = lim [ N-^E^^'>T = [ E{n)T 

N^ooN Jx!^ JviX) 

as desired (using the dominated convergence theorem in the last step, which applies 
thanks to the bound < A + Bd 2 resulting from the Main Assumptions). 

□ 

Next we turn to the asymptotics of the distances, establishing the following key 
property: 

Proposition 2.21. Assume that a sequence vx of symmetric probability measures 
on satisfies 

lim {Sn)*i^n = Si, 

N —^co 

in the distanee topology in W 2 (P 2 (A)). Then any sequenee fj,x such that {6x)*tJ-N 
converges weakly to F S 'P{'P{X)) satisfies 

liminf ^(i(pAr,i/Ar)^ > / d(/r, i/)^r(/r) 

N-ioo N Jv{X) 

and equality holds iff {Sn)*P,n converges to T in the distance topology in W 2 {'P 2 {X)). 
Proof. Consider the isometry 

Sn ■ {X^^\dx(N)) ^ {V{X),dw) {xi, ..,xn) ^ Sn'■= 

defined in terms of the distances. We equip the space V{'P{X)) with the 
Wasserstein (pre-)metric d induced from distance dw onV{X), i.e. we consider 
the subspace IF 2 (P(A)). By Lemma [131 

—d{p,x, vn)"^ = d{{SN)*tJ,N, {Sn)*i^n)‘^ 

We now first assume that {Sn)*P'N converges to F in the d—distance topology in 
W 2 {P 2 {X)). Then the “triangle inequality” for d immediately gives 

lim d{{SN)*fJ-N,iSN)*i^N)^ =d(r,(5y)^. 

N—yoo 

Next we will use the following simple general fact for the Wasserstein distance on 
V{Y,d) : 

d{lJ‘,Syo)'^ = J d{y,yofp.{y) 

which follows from the fact that the only coupling between /i and Sy^ is the product 
p. (Si Syg. Applied to Y = V{X) this gives 

d{{SN)*PN,Syf = f d{p,vYT(p) 

Jv(x) 
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which concludes the proof using that = d(p, u) by the general fact above. 

More generally, if {6N)*^J'N is only assumed to converge to F weakly in 'P{V{X)), 
then the lower semi-continuity of the Wasserstein distance function wrt the weak 
topology instead gives 

liminf 

N^oo N Jv{X) 

Finally, if equality holds above, then, by the previous arguments, 

lim / d{^,v)‘^{5N)^^N = / 

J^^V(X) Jv{X) 

(i.e. the “second moments of {5n)*P^n converge to the second moments of F) and 
then it follows from Proposition 12.21 that {6N)*fJ'N converges to F in the distance 
topology in W 2 {J^{X)). □ 


2.6.2. Conclusion of the proof of Theorem \2.1‘T\. Without loss of generality we may 
set /3 = 1. We start by observing that for any fixed p in V{X) we have, by the 
defining property of , that 

where the rhs converges, by the propositions above, to Jj+i(/r) as N ^ oo, where 
Jj+i{p.) = -I- F{p). In particular, taking p = p-j+i gives 

(2.19) lim sup ) /N < Jj+i {fij+i) 

N—^oo 


where pLj+i is the unique minimizer of Jj+i- 

Next we consider the lower bound. By the minimizing property of we have 
a uniform control on the ^ 2 —distance: 

(2.20) < C 

Indeed, the minimizing property together with the previous bound gives 

<C- 1fW(pWj 

Hence, it is enough to verify that the uniform coercivity property 12.121 holds. But 
this follows the uniform Lipschitz assumption on . 

Now, it follows from the induction assumption and the triangle inequality for 
d that satisfies the assumptions of Lemma 12.21 Accordingly, we may, after 

passing to a subsequence, assume that (J-n ■= dif+i converges as in Lemma r2.18l or 
more precisely that 

iSN)*di^Ji F 

in IFi)^)^)), for some F e W 2 (V{X)). It then follows from Propositions l2. 191 12.201 
and 12.211 that 

(2.21) ImHnf > j dF(p) J.+i (p) 

Combining the previous lower bound with the upper bound 12.191 and using that 
p.j+1 is the unique minimizer of Jj+i then forces F = and 


( 2 . 22 ) 




N—>oo 
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But this means that 


weakly in V{X) and by the equality 12.221 that 

But then it follows from Proposition 12.211 (applied to = <5^^.) that {SN)*^J'N con¬ 
verges to r in the distance topology in W 2 {V{X)), as desired. 

2.7. Convergence in the non-discrete setting: proof of Theorem [T7ll We 

first assume that A > 0. By the Main Assumptions the limiting free energy func¬ 
tional F{fi) := E{p) -|- Hip) is also Lipschitz continuous and A—convex along gen¬ 
eralized geodesics in 'P2IX). Indeed, bv 12.81 Eia) is the limit of the mean energy 
functionals /x 1 —>■ /N which are A—convex along generalized geodesics, 

since E^^'> /N is assumed A—convex (see Lemma EdS]). In particular, by Theorem 
12.141 the gradient flow nt of E emanating from a given /xq S E{X) exists and is 
uniquely determined in the sense of Theorem 12.91 We let Ft := 5^^ be the corre¬ 
sponding flow on 7^2(^ 2 (A)). 

Consider the fixed time interval [0,T] and fix a small time step r > 0. Denote 
by /x'^(t) the discretized minimizing movement of A(/x) with time step r and set 
rj" := 5^t. For any fixed t €]0,T[ we then have, by the triangle inequality, 

d(r^(t),r(t)) < d(rjv(t),r]V(t)) + dim, r"(t)) + d(r]v(0,r"(t)) 

By the isometry property in Lemma 12.31 and the assumed convexity properties 
we have, by Theorem 12.151 that dlT]\flt),T^lt)) < Ct (uniformly in N) and 
(i(r(t), r'^(t)) < tC. Moreover, by Theorem 12.171 limv^.^ dirCfti. T'^D')') = 0 for 
any fixed r. Hence, letting first N ^ 00 and then r —>■ 0 gives limAr->.oo d(rv(t), r(t)) 
0 , which concludes the proof. 

In the case when A < 0 the previous argument still applies (with the error 0 (r) 
replaced by ©(r^/^) according to Remark 12.161 

2.8. Proof of Theorem 11.21 (the static case). By Gibbs variational principle 
(which follows immediately from Jensen’s lemma) is a minimizer of the mean 
free energy E^^\ In particular, for any fixed fj, £ P(]R"') 

(2.23) ^E(^)im) < +1 < C, 

where the last inequality is obtained by taking /x to be any measure with compact 
support. Next observe that by the properness assumption on this gives 

f < AlFW(AxW)-i?<AC-H 

7(Rn)iV A ^ A 

However, on order to apply Proposition 12.201 we would rather need a bound on the 
p—moments for some p > I : 


(2.24) 
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But using the convexity assumption this follows automatically from the bound on 
the first moments using the following well-consequence of Borell’s lemma m (see 
[5^ Appendix III]), which gives a “Reversed Holder’s inequality”: 

Lemma 2.22. Let ^ be a log concave measure, i.e. p = e~'^dx for some convex 
functions cj). Then, for any q > p : 



where the constant Cp^q^n is independent of p 

To prove the bound 12.241 we first observe that for any p 





N 


Jr^ 


where denotes the “first marginal” of i.e. its push-forward under the 

natural projection (R")^ —>• R.” onto the first factor. But, by assumption 
is log concave and hence, by the Prekopa inequality [5^, so is its first marginal 
Applying the previous lemma thus gives 


JR" 

for any p > 1 where Cp^n is independent of N. In particular, by Proposition 12.21 
combined with the isometric embedding in Lemma 12.31 we get that '■= 

converges to some P in the distance topology in Wp{V{X)) for any p £ [1, cx)[. 
Applying this to p = I and invoking Proposition [2T9] thus gives the convergence of 
the mean energy: 

hm 1 [ Eip)T 

AT-j-oo iV JxN J'P(X) 

We then deduce, using the asymptotics of the entropy in Prop 12 .1^ precisely as in 
the proof of Theorem 12.171 that 


[ F(p)r(p) < hminf IfW(pW) < 

J-p{X) Af-i-oo N N 

Next we get, using the upper bound l^T^ and Proposition l2.19l applied to px = , 

that 

[ E{p)r < hminf IfW(pW) = F(/i*) = inf R(p) 

J'P(X) Af—>00 A W2iX) 

Hence, it must be that P = (5^, which concludes the proof of the convergence 
assuming the existence and uniqueness of the minimizer p,. In fact, the existence of 
p, also follows from the previous argument: indeed, since, by well-known properties 
of the entropy |2H], we have E[^^\p)/N > E[{p[^'^) the previous argument gives 
that any limit point p* of the sequence in 'P{X) (which exists by tightness 
and, as explained above, is in Wp(X) for any p) satisfies 


F(p.) < hminf IfW( pW) = F(p,) = inf R(p) 
N-¥oo iV W2(X) 
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and hence minimizes F on W 2 {X). Finally, since fj,\ is log concave so is the 
limit As for the convergence 11.101 it follows immediately from the formula 
— ^logZ^ = and the proof of the above does not use neither exis¬ 

tence nor uniqueness of a minimizer of F. 


3. Permanental processes and toric Kahler-Einstein metrics 


3.1. Permanental processes: setup. Let P be a convex body in K" containing 
zero in its interior and denote by pp the corresponding uniform probability measure 
on P, i.e. P — lpdX/V{P), where dX denotes Lebesgue measure and V{P) is the 
Euclidean volume of P. Setting Pk := Pr\{7j/k)'^, we let Nk be the number of points 
in Pk and fix an auxiliary ordering of the Nk elements of Pk- Given a 

configuration {xi, of points on A := R" we set 

(3.1) := ilog V 

rC 

which, as explained in the introduction of the paper, can be written as the scaled 
logarithm of a permanent. To simplify the notation we will often drop the subscript 
k and simply write Nk = A, since anyway A ^ oo iff fc —>■ oo. We will denote by 
v) the Monge-Kantorovich optimal cost for transport between the probability 
measures p and with respect to the standard symmetric quadratic cost function 
c{x,p) = —X ■ p : 


(3.2) 


C(p, p) := inf — [ 
'I iA 


x-jry, 


XxX 


where the 7 ranges over all couplings (transport plans) between p and v. 


Proposition 3.1. The Main Assumptions for are satisfied with A = 0 and 

E(fi) = —C{p, vp)i Equivalently, formulated in terms of the Wasserstein L'^ —distance 


(3.3) 




/ 


+ Cp, cp := 


1 


/ 


P^Vp 


In particular, —C{-,vp) is convex along generalized geodesics. 


Proof. This follows from the results in [S]. In fact, the first and second point follows 
immediately from basic fact that if fig is a family of smooth convex functions on 
K" and v a probability measure on the parameter space, then (j) := log/ e^’‘v{s) 
is also convex and V(/) is contained in the convex support of {V^s}, which in the 
present case is contained in kP. Hence, € P which is uniformly bounded, 

since P is a convex body and in particular bounded. Finally, the convergence of 
E(n) 

was shown in [S] for pL with compact support (which is enough by Lemma 
12.51) . The convexity of —C(-, vp) then follows from Lemma [2.131 Equivalently, this 
means that —Xdw 2 {hj ppY is —1—convex. In fact, as shown in [2] using a different 
argument is —1—convex for any fixed p G P 2 (R")- GI 


Next, we recall that the Monge-Ampere measure MAffi) of a convex function 
on R." is defined by the property that, for a given Borel set E, 


[ MA{4>) := [ 
Je J(i 


dX, 


{d4’KE) 
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where dX denotes Lebesgue measure and d4> denotes the subgradient of (f) (which 
defines a multivalued map from M" to K"). In particular, if S then 

= det{d^4’)dx, 

where d^(j) denotes the Hessian matrix of (j). We will denote by Cp the space of all 
convex functions (f) on M" whose subgradient d(j) satisfies 

C P 

and we will say that (j) is normalized if (/>(0) = 0. By the convexity of (/) the gradient 
condition above equivalently means that (j) grows as most as the support function 
of P, where (/'p(a;) := suppgpp • x 

By Brenier’s theorem m, given fj, = pdx in P 2 (R"') there exists a unique nor¬ 
malized (f> S Cp such that 

(3.4) MA{4>) = pdx, 

which equivalently means that the corresponding map Vcj) from R." to P sat¬ 
isfies 

{V4>)*p = vp 

Given the previous proposition we can use the differentiability result in [2] for the 
Wasserstein distance to get the following 

Lemma 3 . 2 . The minimal subdifferential of —C{-,vp) on the subspace 7 ^ 2 ,ac(®.") 
of all probability measures in 7^2 (R”) which are absolutely continuous wrt dx, may, 
at a given point pdx, be represented by the L°° —vector field where (j) is the 

unique normalized solution in Cp to the equation \3.4\ 

Proof. Given formula [3.3l this follows immediately from Theorem 10.4.12 in [2] and 
the fact that if 7 ^ 2 ,ac(R")j then Brenier’s theorem gives that the optimal transport 
plan (coupling) from R" to P realizing the infimum defining dw 2 (/^j ^p)"^ is given by 
the map V<f>, where 4> solves the eauation l3.4l Since the barycentric projection 
appearing in Theorem 10.4.12 in |2] for the transport plan defined by a transport 
map gives back the transport map (see [21 Th, 12.4.4] this concludes the proof. See 
also |5] for a direct variational proof of Brenier’s theorem which can be seen as the 
real analogue of the variational approach to complex Monge-Ampere equations in 

a. □ 

3.2. Existence of the gradient flow for Fp{p,). Given (3 g]0, oo] we set Fp{p) := 
-C{p,vp) + H{p)/j3 

Proposition 3 . 3 . The gradient flow pt of Fp on 7 ^ 2 (R") emanating from a given 
po exists for any (3 g]0,oo]. Moreover, for f3 < oo we have that pt = Pt{x)dx 
where pt has finite Boltzmann entropy and Fisher information and p{x,t) := pt{x) 
satisfies the following equation in the sense of distributions on ]R"x]0,oo[ 

(3.5) ^ = ^Apt + V • {ptVcft), 

where (ft is the unique normalized solution inCp to the equation \8.4\ and 'S/eft defines 
a vector field with coefficients in 

Proof. Given the previous lemma this follows immediately from Thm 8.3.1 and Gor 
11.1.8 in |2] (the case (3 = oo has previously been considered by Brenier [201 HU [H] 
by lifting the problem to the space of maps from R." to R" where Hilbert space 
techniques can be applied ). □ 
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More generally, as explained in the introduction it is natural to introduce a 
parameter 7 G [ 0 , 1 ] and a back-ground potential V{x)^ i.e. a convex function on 
R" satisfying the growth condition 1 1.1 81 Then one replaces with its weighted 
generalization E^y defined by formula 11.191 Then the previous proposition still 
holds with Fp replaced by the corresponding functional and (j)t in the evolution 
equation 13.51 is replaced by 7(^4 -|- (1 — 7 )!^ and with /3 = 1 (up to rescaling time t 
and the potential V this is equivalent to taking /I = 7 ) 

3.3. The dynamic setting: Proof of Theorem 11.31 By Proposition 13.11 the 
Main Assumptions are satisfied. 

3.4. The static setting: Proof of Theorem 11.41 Let us first verify the proper¬ 
ness assumption in Theorem 11.21 holds for 7 < Rp. To this end we first assume that 
bp = 0 and observe that 

N 

(3.6) ....,xn) > - 0 ( 1 )^ \xi \ - 0 ( 1 ) 

Indeed, applying Jensen’s inequality to the concave function log gives 

= ylog^ V - I log ^ > 

k Nk'. ^ k Nk'. 

1 ^ 

i=i O’ 

But for any fixed i we have that J2aP<y{-i) ~ ~ ^)-SpjGPnz/fePi hence we 

get a Riemann sum: 

^ Pj :=bP =bp + o{l), 

where bp := fppup, which is assumed to vanish and hence the inequality 13.61 
follows. But then the properness for E^^y (defined by formula 11.191) in the case 
bp = 0 follows immediately from the definition of E^^ and the growth assumption 
ll.l 8 l on V ensuring that that V(x) > \x\/C — C since 0 is assumed to be an interior 
point of P. Finally, the case then bp ^ Q can be reduced to the previous case by 
translating P. More precisely, by the previous argument 

N 

+ 0 ( 1 ) > ^ (7 + o{l))x, ■bp + {l- 7 )R {xij) 

1=1 

But, as shown in the proof of Theorem 2.18 in [5] Rp (defined by formnla 11.211) is 
the sup of all r £ [ 0 , 1 ] such that rx ■ bp + {1 — r)(j)p{x) > 0 and since |(/'p(a:) — 
V{x)\ < C and 7 < Rp this gives the desired properness. The convergence of the 
Boltzmann-Gibbs measures, as N -^ 00 , now follows from Theorem 1 1.2 1 Moreover, 
the convergence as t ^ 00 for IV fixed follows from well-known resnlts about the 
linear Fokker-Planck equation with a convex potential E such that J e~^dx < 00 
(see for example [14] and reference therein). 

Finally, to prove the divergence of the partition function for 7 = we first recall 
that if '0 is a convex function on then a necessary condition for J e~'^dx < 00 
is that ^(a;) ^ 00 as \x\ —>■ oo(as follows, for example from Borell’s lemma |17l 
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Lemma 3.1] which gives that the integrability also holds for ■0 — ejxj, for e any 
sufficiently small number ini Theorem 3.1], and hence, by Jensen’s inequality, 
^(x) — ejxj > —Ce := — log / To violate the previous condition it is 

clearly enough to find a vector a G such that 1 1 —>■ '!/i’(ta) is an affine function on 
K. We now consider the convex function tp := on observe that 

for any fixed a G R” we can write, with a := (a, a,...., a), 

(3.8) i/;(a)-'!/)(0) = ((1 - i?p)(^p(a) + i?pa ■ 5p) + a-i5fc, 4 := i?p(6p^-6p) G M" 

(compare formula 13.71) . Moreover, as shown in the proof of Theorem 2.18 in [5] 
when a is taken as a normal vector to a facet of P containing the point q appearing 
in formula [1.2 II) the bracket in formula [3^ vanishes. But this means that 1 1 —>■ 
is an affine function on R and hence f e~'^dx = oo. Since \(j)p{x) — V{x)\ < C the 
divergence also holds when (pp is replaced by V, which concludes the proof. 

3.5. The tropical setting. The results above are also valid when the permanental 
interaction energy E^^'‘'>{xi, ...,xn^) is replaced by its tropical analog, i.e. the 
following convex piecewise affine convex function 

EtZp{Xl,-,XNk) ^ - \-XN-Pa{Nt,) 

In other words this is a tropical permanent, i.e. the permanent of the rank N matrix 
(xi ■ pj) in the tropical semi-ring over R (i.e. the set RU {—oo} where the plus and 
multiplication operations are defined by maxja, 6 } and a -I- 6 , respectively [44| 1. 
Equivalently, in terms of discrete transport theory this means that 

Eirop{xi,...,XN^) := -C{{5N{x),5Nip))) 

Passing to the tropical setting has, in particular, computational advantages. Indeed, 
while all known methods for evaluating (general) permanents take exponential time, 
the tropical permanent above is, by its very definition, the optimal value of a 
linear assignment problem and can be computed using an algorithm of cubic time 
complexity (see the discussion in [24] 1. 

4. Outlook 

In this final section we point out some relations between the limiting evolu¬ 
tion equation appearing in Theorem [L3] (whose static solutions correspond to toric 
Kahler-Einstein metrics) and other well-known evolution equations. We also indi¬ 
cate some relations to sticky particle systems appearing at the microscopic level (i.e. 
for finite N) and the complex geometric picture. These relations will be elaborated 
on in a sequel to the present paper HU. 

4.1. Relation to other evolution equations and traveling waves. In the one¬ 
dimensional case when P := ]—a_, —a+] integrating the evolution equation for pt 
in Theorem [T3| once reveals that the bounded decreasing function u(x, t) := —dx4>t 
(physically playing the role of a velocity field) satisfies Burger’s equation |43| with 
positive viscosity k := j5~^ : 

dtu = nd^u — udxU 
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with the left and right space asymptotics limj^-^j-oo w(x, t) = a±. We recall that 
Burger’s equation is the prototype of a non-linear wave equation and a scalar con¬ 
servation law (which is used, among many other things, as a toy model for tur¬ 
bulence in the Navier-Stokes equations [S^). Interestingly, the barycenter bp of 
the polytope P coincides, in this one-dimensional situation, with minus the speed 
s := (a+ + fl -)/2 of the time-dependent solution u in the terminology of scalar 
conservation laws [50]. Hence, the vanishing condition bp = 0, which in general is 
tantamount to the existence of a stationary solution pt = p (as discussed in con¬ 
nection to Theorem [E3D simply means, from the point of view of non-linear wave 
theory, that the speed s vanishes. 

Similarly, the function (j){x,t) := which in complex geometric terms is 

a Kahler potential, satisfies (after the appropriate normalization) the following 
viscous Hamilton-Jacobi equation (known as the deterministic KPZ-equation in 
the literature on growth of random surfaces mmo]): 

(4.1) dtcj) = Kdlcj) + . 

In the general higher dimensional case the evolution equation 11.141 (which is dif¬ 
ferent than the higher dimensional version of Burger’s equation) can be seen as 
a dissipative viscous/diffusive version of the semi-geostrophic equation appearing 
in dynamic meteorology (see [521 [1] and references therein and | 2 T| for a similar 
situation in cosmology). Moreover, since 

EifJ-) = ’^p) + \ J + C, 

where d denotes the Wasserstein distance, the evolution equation II. 141 can also 
be seen as a quadratic perturbation (with diffusion) of the “geodesic flow” on the 
Wasserstein space (compare |2] Example 11.2.10]), which in the one dimensional 
case appears in connection to the Sticky Particle System m- As will be shown 
in [12] the large time asymptotics of the fully non-linear evolution equation 11.141 
for the probability density pt in R" are governed by traveling wave solutions in R." 
whose speed coincide with minus the barycenter bp oi the convex body P : 

Pt{x) = p{x — bpt) + o{t), t —>■ 00 

where the error terms o{t) tends to zero in L^(R") (and even in relative entropy) and 
where the limiting profile p is uniquely determined from a variant of the Monge- 
Ampere equation 11.171 together with the condition that its barycenter coincides 
with the barycenter of the initial data (thus breaking the translation symmetry). 
In complex geometric terms p corresponds to a certain canonical Kahler-Einstein 
metric w on AT with conical singularities “at infinity”, playing the role of Calabi’s 
extremal metrics in this context. More generally, as will be elaborated on in [12], the 
results above apply in a more general setting where the measure z/p is multiplied 
by a density g, which amount to replacing the Monge-Ampere equation 
with g{'S/(j))MA{(j)) and which from the point of view of scalar conservation laws 
corresponds to a general concave flux function / (when n = 1 ). 

4.2. The microscopic picture: sticky particles in R". It can be shown that 
the attractive Newtonian interaction energy in R is the one-dimensional version of 
the tropical permanental energy ...xn) appearing in Section 13.51 In the 

general higher dimensional setting it turns out that a very concrete interpretation 
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of the corresponding EVI gradient flow of ...xn) on can be given; in 

particular the particles perform zigzag paths with velocity vectors contained in the 
polytope —P generalizing the sticky TV— particle system on the real line. ^Moreover, 
there is a static solution to the corresponding deterministic TV— particle system if 
and only if the “discrete” barycenter of P vanishes: 

^{Pi + ■■■+Pn) = 0 

which is consistent with the fact that the discrete barycenter can be interpreted as 
the mean velocity of the particles. In general, any initial configuration of points 
(cci, ...a:Ar)(0) is assembled, in a finite time, into a single particle a;*, namely the 
barycenter of {a;i, ...xn} which moves at the mean velocity above. The results in 
the present paper can also be used to study the large TV— limit of this deterministic 
system (which can be seen as a dissipative version of the Hamiltonian particle 
system introduced in [53] as a discretization of the semigeostrophic equations). 
But the key point of our approach is that it allows noise to be added to the particle 
system. Then the role of the large TV— limit of a;* is played by the volume form 
/r* of a Kahler-Einstein metric on the toric variety determined by the polytope P 
(compare the discussion in Section HTTl) . 

Interestingly, a similar particle system on M" appears in Brenier’s approach m 
122 ) to the early universe reconstruction problem in cosmology [36| (in connection to 
the so called Zeldovich approximation). In fact, our results can be used to validate 
the formal large TV— limit of the TV— particle system with noise introduced in [221 
Section 2.3]0Details will appear in m- 

4.3. The complex geometric picture. In this final section we provide some 
complex geometric motivation for the present paper - a more detailed account, 
including the relations to the Yau-Tian-Donaldson conjecture, will appear elsewhere 

Let X be an n—dimensional compact complex manifold. A metric g on X is 
said to be Kahler-Einstein if g has constant Ricci curvature and g is Kahler, i.e. in 
local holomorphic coordinates g can be represented as the real part of the positive 
definite complex Hessian of a local function (j){z) called the Kahler potential 

of g. If such a metric g exists with positive Ricci curvature, then X is necessarily 
a projective algebraic variety which is Fano, i.e. the holomorphic (anti-canonical) 
line bundle L := det(TA) over X is positive. 

As shown in [10] a Fano manifold comes with a sequence of canonical TV— particle 
random point process. The number of particles TV arise as the pluri-antigenera of 
A : 

N = Nk :=dimH°(A,L®'=), /fc = 1, 2,3,.... 

where iL°(A, L®^) denotes the complex vector space consisting of the global holo- 
morphic sections of the k th tensor power of L. The Fano condition ensures that 
Nk —>■ oo as A: ^ cx). The local density of the corresponding canonical symmetric 


^When n = 1 the dynamics is determined by the property that total mass and momentum 
is conserved in collisions and that the particles stick together when they collide; see m and 
references therein. 

^ As pointed out in |22l Section 2.3] the formal argument used there, which is based on the 
classical Freidlin-Wentzel theory (as in 1311 ). would require a Lipschitz bound on the drift (while 
only a one-sided bound is available). 
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probability measure on X^’’ is defined by 

(4.2) 

:= \det{zi,...ZN,)\~‘^^^ , det(zi, ...zjvJ := det(s*(zj)), 

where det(zi, , L®^) is the Vandermonde type determinant formed 

from a given base si,...sisi^ in i/°(V, L®^) and is the corresponding normal¬ 
ization constant ensuring that the probability measure has unit mass (by homogen- 
ity is independent on the choice of base). However, since the local density 

, ZNk ) has singularities (for example when two points on X merge) the 
normalization constant may be infinite, which means that the random point 
processes are only well-defined if < oo. Such a Fano manifold X was called 
Gibbs stable in [10], where it was shown that the condition can be rephrased in 
purely algebro-geometric terms (see also [37] for further developments). It was 
conjectured in uni that this condition is equivalent to X admitting a (unique) 
Kahler-Einstein metric (which necessarily has positive Ricci curvature) whose vol¬ 
ume form may be recovered as the deterministic large V—limit of the empirical 
measures of the corresponding random point processes. 0 

The motivation for the present paper comes from a dynamic approach to the 
latter conjecture where one introduces the interaction energy 

E'^^'^'>{zi,...,zn) := ■ilog|det(2;i,...2;jvJP 

which is attractive, in the sense that it tends to —c» as two particles merge. Locally, 
this object is represented by a plurisubharmonic function, but in order to get a 
globally well-defined function on X^’^ one also has to fix a back-ground Kahler 
metric g on V (representing the first Chern class of X) whose volume form dVg 
then induces a metric H-H on L which is used to replace the absolute values above. 
The point is that, if X is Gibbs stable, the canonical probability measure 
on XNk 

can then be represented globally as the corresponding Gibbs measure at 
inverse temperature /3 = 1 (which is independent of the choice of metric g) : 

V ZNk " / 

i.e. as a determinantal point process on X at negative temperature. The different 
zero-temperature case was studied in j^. 

At any rate, even if X is not Gibbs stable one can still look at the stochas¬ 
tic gradient flow of on the fold product of the Riemannian manifold 

{X,g). From this dynamic perspective Gibbs stability simply means that the cor¬ 
responding stochastic process has an invariant measure, to wit, Accordingly, 

the natural dynamic generalization of the conjecture referred to above is that a 
(unique ) Kahler-Einstein metric gK e exists precisely when the stochastic gradient 
flow of £(^'=1 admits a stationary measure and then its volume form can 

be recovered from the joint large N and large t limit of the flow. More precisely, 
conjecturally the large V—limit of the corresponding stochastic gradient flows is 


^The convergence of the processes in the opposite case when the dual det(T*X) of det(rX) 
is positive was settled in HQ] (the limit is then the volume form of the unique Kahler-Einstein 
metric on X with negative Ricci curvature, whose existence was first established in the seminal 
works of Aubin and Yau). 
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described by the complex version of the evolution equation [331 obtained by replac¬ 
ing the real Monge-Ampere operator with its complex counterpart. The latter flow 
is, at least formally, the Wasserstein gradient flow of a free energy type functional 
F{^) on 7^2 (A", g) and F can be identified with the K-energy functional on the space 
of Kahler metrics in ci(X) (using the Calabi-Yau isomorphism) [7]. Unfortunately, 
the study of the latter flows is plagued by various analytical difficulties stemming 
from the singularities of and the lack of convexity. For example, even in the 

simplest case when X is the Riemann sphere, i.e. the one-point compactification 
of the complex plane C, so that is simply the attractive logarithmic pair 

interaction between Nk equal charges on C, the convergence of the large Y—limit, 
for a fixed time, is a long-standing open problem (however, see |35| for very recent 
partial results). 

4.3.1. The toric setting and its tropicalization. The complex geometric setting which 
is relevant to the present paper appears when A is a toric Fano manifold, i.e. X 
admits a holomorphic action of the the real n—torus T such that (A, T) can be real¬ 
ized as an equivariant compactification of the complex torus C*” (with is standard 
T—action) [34]. Such a compactification A is determined by a convex polytope P, 
which has the property that under the dense embedding of C*" into A the com¬ 
plex vector space £°(A, L®^) may be identified with the space of all holomorphic 
Laurent polynomials f{z) on C*" of the form 

f{z) = Y. 

(using multindex notation). In particular, introducing an ordering TOi,...mjVj, on 
the integer points of kP fl Z" gives a bases Smiiz),SmM^. of multinomials in 
H'^{X, which can be used to represent 

(4.3) = y: 

ctGSn 

Now, the real vector space MF makes it appearance when introducing logarithmic 
coordinates on C*", i.e. as the image of the Log map 

Log: z^a;:=(log|zi|2,...,log|z„n, 

whose fibers are the orbits of the action action of T. Using this map T—invariant 
metrics on L —>■ A with positive curvature may be identified with convex func¬ 
tions (j){x) on K" such that (9(/))(K”) C P. In this picture the permanental den¬ 
sity Per(a;i,..., xat^) arises as the push-forward to K", under the Log map, of the 
determinant density 14.31 In other words, the smooth convex permanental energy 
Ep^r (xi, ...xn) (formula ll.121) on R" is an averaged version of the singular plurisub- 
harmonic interaction energy on C*" : 

(4.4) E^^Xxu -Xn) = Uog[ 

Similarly, its tropical version E^^^pixi, ...xn) is the piecewise affine convex function 
on R”^ obtained as the tropicalization of the Laurent polynomial det(zi, ...zn^) on 
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^*nNk a Accordingly, Theorem 11.31 and Theorem 11.41 should be seen in the light of 
the well-known philosophy of replacing an elusive complex geometric problem by 
a more tractable convex geometric one, by the process of tropicalization (see, for 
example, m). 

From the complex geometric point of view the equation 11.201 is precisely the 
one appearing in Aubin’s continuity method for Kahler-Einstein metrics. In the 
complex setting the sup over all 7 S [ 0 , 1 ] for which the corresponding complex 
Monge-Ampere equation has a solution coincides with the the differential geometric 
invariant R{X) of X defined as the greatest lower bound on the Ricci curvature of 
any Kahler metric uj € ci{X) [^. Moreover, as first shown in [51] the invariant 
Rp of P coincides with R{X) when X is the toric Fano variety corresponding to 
P. In this complex geometric context the last point in Theorem 11.41 implies one of 
the inequalities in the conjectural equality relating R{X) to the algebro-geometric 
invariant 7 (A) introduced in |10| . when A is a toric variety. To explain this we 
recall that 


(4.5) 7 (A) := liminf lct(X>fc/fc) > 1, 

k—¥oc 

where \ct{'Dk/k) is the log canonical thresholds of the divisor Vk in X^’‘ cut out 
by the holomorphic section of corresponding to det( 2 i, ...zn^)- The number 

lct(I?fc/fc) may be analytically defined as the sup over all positive numbers 7 such 
that ||det||~^^'^^ is hnite. 


Corollary 4.1. Let X be a toric Fano variety. For any positive integer k 

^lct{Vk) < R{X) 

In particular, on any toric Fano variety X the following inequality holds: 

i?(A) > min{ 7 (A), 1} 


Proof. Write the complex coordinates on C*" as where G R 

and y°‘ G ]R/[0, 27r]Z for a = l,..,n. In particular, x = logz in the notation of 
Section 14.3.11 Representing dV = e~^ dx A dy for a T—invariant function V 

on C*", which we identify with a function V{x) on R" such that (cl(/))(R") C P, 
gives 


/ if 

\Jy 


1 


(|det(zi,...,2:ivjp) "(.^ ^ 

yG[0,27r["-^fc 




where H-H denotes the metric on induced by dV. Hence, applying Jensen’s 

inequality to the convex function /(s) = on ]0, oo[ and using the relation [4.41 

gives 

[ \\det{zi,...ZNk)\n^^''dV®^>‘ > [ Per (xi,..., xatJ”* 

J 


^Incidentally, tropicalization may be interpreted as a zero-temperature limit by writing the 
tropical sum max{a, b} as the limit of T~^ log(eT“ + as T ^ 0; compare the discussion in 
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But, according to Theorem 11.41 the latter integral diverges if 7 > Rp and that 
concludes the proof. □ 


The last statement in the previous corollary can be viewed as a toric complement 
to Fujita’s very recent result m, which in particular gives that if ^{X) > 1 on a 
Fano variety X, then X is K-semistable and hence R{X) > 1 (which in the toric 
case is equivalent to the existence of a Kahler-Einstein metric mi This gives an 
illustration of how the process of tropicalization can be used to gain information 
on the original complex geometric situation. 


5. Appendix: The Otto calculus 

In this appendix we briefly recall Otto’s |62| beautiful (formal) Riemannian in¬ 
terpretation of the Wasserstein metric c ?2 on 7^^(K”). The material is included 
with the non-expert in mind as a motivation for the material on gradient flows on 
recalled in Section [231 


5.1. The Otto metric. For simplicity we will consider probability measures of 
the form p = pdx where p is smooth positive everywhere (in order to make the 
arguments below rigorous one should also specify the rate of decay of p at c» in 
K"). The corresponding subspace of probability measures in P 2 (K.") will be denoted 
by V. First recall that the ordinary “affine tangent vector” of a curve pt in V at 
p := po, when p* is viewed as a curve in the affine space L^(R"')) is the function p 
on R” defined by 


p{x) 


dptjx) 
dt |t=o 


Next, let us show how to identify p with a vector field Vp in ^^(pdx,R"), which, by 
definition, is the (non-afhne) “tangent vector” of pt at p, i.e. Vp € TpV. First, since 
the total mass of pt is preserved we have / pdx = 0 and hence there is a vector field 
V on R" solving the following continuity equation: 


(5.1) 


p = -V ■ (pv) 


In geometric terms this means that 

(5.2) ptdx = {F^)^{podx) + o{t), 

where is the family of maps defined by the flow of V. Now, under suitable 
regularity assumptions Vp may be defined as the “optimal” vector field v solving the 
previous equation, in the sense that it minimizes the norm in L^{pdx, R"). The 
Otto metric is then defined by 


(5.3) 


g{vp,Vp) = inf J p\v\‘^dx = j p\vp\^dx, 


which can be seen as the linearized version of the defining formula 12.61 for the 
Wasserstein metric. By duality the optimal vector field Vp may be written as 
Vp = for a unique normalized function (j) on R" (under suitable assumptions). 
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5.2. The microscopic point of view. Let us remark that a simple heuristic 
“microscopic” derivation of the Otto metric can be given using the isometry defined 
by the empirical measure Sn (Lemma 12.31) . Indeed, given a curve (xi(t), ...^XN{t)) 
in the Riemannian product with tangent vector 

t = 0 we can write its squared Riemannian norm at (a;i(0), ...,a:jv(0)) as 

2 „ 


(5.4) 


dxi{t) dxi(t) 


dt 


dt 




where 5M{t) '■= and v is any vector field on X — R" such that v{xi) = 

|t-o' setting pt := (5iv(t) the vector field v satisfies the push-forward 

relation 15.21 (with vanishing error term). Moreover, since passing to the quotient 
/Sn does not effect the corresponding curve pt, minimizing with respect to the 
action of the permutation group Sn in formula 15.41 corresponds to the infimum 
defining the Otto metric in formula 15.31 


5.3. Relation to gradient flows and drift-diffnsion equations. If G is a 

smooth functional on V then a direct computations reveals that its (formal) gradi¬ 
ent wrt the Otto metric at p corresponds to the vector field v{x) = In 

other words, the gradient flow of G{p) may be written as 

(5.5) - - = Va; • {pvt[x)), Vt{x) = Va;—-- 

dp 

In particular, for the Boltzmann entropy H{p) (formula l2.1l) one gets, since = 

logp (using that the mass is preserved), that the corresponding gradient flow is the 
heat (diffusion) equation and the gradient flow structure then implies that id(pt) 
is decreasing along the heat equation. Moreover, a direct calculation reveals that 
H is convex on V in sense that the Hessian of H is non-negative and hence it also 
follows from general principles that the squared Riemannian norm \XH\^{pt) is 
decreasing. In fact, by definition |V7Lp(p) coincides with the Fisher information 
functional I{p) (formula 12.1|) . More generally, the gradient flow of the Gibbs free 
energy Pp is given by the diffusion equation with linear drift VxV : 

(5.6) % = ^AxPt + V, • {ptVxV), 

at p 

often called the linear Fokker-Planck equation in the mathematical physics litera¬ 
ture. The study of the previous flow using a variational discretization scheme on 
P^(R") was introduced in [IS] (compare Section |S|). 
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